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ABSTRACT 


A generalized automated data acquisition system was designed for the Naval Postgraduate 


Sct ol Low Speed Wind Tunnel. A specific application of this system was to upgrade the current 


"Cylinder Drag Experiment" conducted during AA2801 “Aero Laboratories I“, an introductory 


aeronautical laboratory course taught at the Naval Postgraduate School. Two methods of drag 


determination were used: pressure distribution and wake analysis (momentum method). Data from 


these two methods were collected by a system based on a high speed analog/digital computer board, 


a standard 486 IBM-type PC and data acquisition software. Characteristic methods of reducing data 


from this experiment are discussed. The results obtained by analyzing the acquired data compared 


favorably to empirical data from previous circular cylinder coefficient of drag experiments. This 


automated data acquisition system will facilitate future research and instructional use of the wind 


tunnel. 
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TABLE OF SYMBOLS AND ABBREVIATIONS 


total (stagnation) pressure 

Static pressure 

dynamic pressure = pp - p 

freestream density 

freestream total pressure (associated with freestream velocity) 

freestream static pressure (associated with freestream velocity) 

freestream dynamic pressure = (Po)o - Po 

wake total pressure measured behind the cylinder at the probe 
((Po)y assumed to = (py). when taken outside wake) 

static pressure measured behind the cylinder at the probe 
(p,, assumed to = p,, when taken outside wake) 

wake dynamic pressure = (Po), - Pw 

pressure measured at the port on the cylinder 

pressure coefficient 

drag coefficient (3-D) 

drag coefficient relating to total drag or profile drag 

drag coefficient -elating to form drag only 

uncorrected drag coefficient 

azimuth of the rotator (signal from potentiometer) 

position of the traverse (signal from potentiometer) 

radius of cylinder 

gas constant 

diameter of the cylinder = 2x R 

incremental area 

incremental pressure drag force 

incremental form drag component 

incremental induced drag component 

drag (2-D) 

length 

width 

height 

cross-sectional area 

measured force 

reference chord 

reference area 

area (2-D) 

volume (3-D) 

blockage factor 

freestream velocity 

velocity as measured along the traverse behind the cylinder 


the velocity into a closed system 

the velocity out of a closed system 

total drag unless specifically defined as “profile drag” 
shear drag - drag due to skin friction 

form drag - drag due to separation 

Reynolds number 

(subscript) effective (corrected value) 

d (subscript) uncorrected value 


hu wi tb th one 


Conversion Constants 


latitude correction = -0.0245 inches Hg 

temperature correction (see Table 4-2 on bulletin board near tunnel) 
Hai = 3.8 x 107 (Ibf-s)/ft 

R,,, = 53.3 (ft-lbf)/(lbm-°R) 

°R = °F + 460 

1 ft = 30.48 cm 

1 cm H,O = 0.0142 psi 

1 psi = 2.036 inches Hg 

Ghest section = (Ap - 0.243)/.895 (NPS wind tunnel calibration equation) 
TF = 1.04 (NPS wind tunnel turbulence factor) 

g, = 32.174 ft/sec? (universal gravitational constant) 
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I, INTRODUCTION 
A classic experiment, drag measurement over a circular cylinder in a wind tunnel, 


serves as an excellent primer on wind tunnel use, data acquisition and basic 


aerodynamics. This thesis served to upgrade an existing AA2801 (Aero Laboratories I) 


experimental lab (see Figure 1) while further exploring the theory behind the methods 
used to measure drag. Chapter II examines the Naval Postgraduate School’s Low Speed 
Wind Tunnel, including a detailed procedure for calibrating the wind tunnel. Chapter 
II will discuss the theory behind drag measurements around a circular cylinder. Chapter 
IV compares the former and new procedures used in the AA2801 laboratory. The results 
of the new experiment using actual class data are then analyzed in Chapter V and some 
conclusions are drawn in Chapter VI. Appendix A contains the former AA2801 Cylinder 
Drag Experiment Handout. Appendix B contains a copy of the upgraded handout that 
was used on 9 November 1993 with a class of graduate students. Appendix C contains 
a readout of the experimental data obtained. Appendix D is the MATLAB script file 
"DRAGI.m" used to analyze data from the upgraded lab. Appendix E has several 
MATLAB script files for creating velocity tables and curves which were derived from 


wind tunnel calibration data. 


Figure 1. The Test Section Showing The Wake Traverse and Probes Apparatus 


i. N.P.S. LOW SPEED WIND TUNNEL 


A. DESCRIPTION AND OPERATING PROCEDURES 

The NPS Laboratory Manual for Low Speed Wind Tunnel Testing [Ref. 1] 
provides an adequate discussion of the basic system description of the NPS Low Speed 
Wind Tunnel and its chief operating procedures and should be read thoroughly by all 
users. What follows is a summary of this description and those procedures of interest 
to the Cylinder Drag Experiment. 


Figure 2 depicts the layout of the NPS Low Speed Wind Tunnel and the location 


of the test section. The wind tunnel is a closed-circuit, single-return type. It is powered 


by a 100 hp electric motor which drives a three-blade variable pitch fan. 

A four-speed common truck transmission is used in-line with the motor and allows 
the wind tunnel to achieve tunnel speeds up to 200 MPH. The Cylinder Drag 
Experiment in AA2801 operates at approximately 115 MPH and for this, fourth gear is 
used. The contraction cone prior to the test section has a 10:1 contraction ratio which 
serves to increase the velocity of the airflow as well as causing the flow to become more 


uniform entering the test section. 


! Test Section 


64.1 Feet 


Figure 2. NPS Low Speed Wind Tunnel 
The dimensions of the test section are given in Table 1. 


TABLE 1. Test Section Dimensions. 


total cross sectional area (C) 
above the reflection plate: 


height (hi) 
cross sectional area (C.,) 


Figures 3 and 4 are views from inside the test section. The wind tunnel uses a 
water micromanometer to reference tunnel velocity by measuring Ap, the difference 


between two static pressures: p, at the beginning of the contraction cone and p, at the exit 


Figure 3. Test Section View - Into the Freestream 

of the contraction cone. Since the airspeed is increased between the beginning and exit 
of the cone, there is a decrease in the static pressure as the air seeks to returm to the 
lower ambient static pressure. This equalization is due to a breather slot, which allows 
ambient air to enter the circuit and mix with air leaving the test section. Hence, Ap 
changes primarily with increasing (or decreasing) p,, aS p, remains essentially constant. 
Figure 5 shows the water micromanometer. Through calibration and associated equations 


and charts (see Section C), this Ap is related to airflow velocity. 


Figure 4. Test Section View - Looking Down Wind 


Table 2 illustrates the operating instructions used during this lab. 


After confirming micromanometer reads 0.00 with tunnel 
fully stopped, set selected Ap in window 


Check tunnel doors closed and test section free from FOD; 
check transmission for proper gear; turn power switch to ON 
position 
Depress and hold START button until green light goes out 


momentarily (meaning fan is turning and set to minimum 
pitch); release START button 


Wait 5-10 seconds; depress COARSE UP button until water 

rises to enter the crosshair ring on micromanometer; toggle 

FINE UP/DOWN to maneuver meniscus to rest on crosshair; 
tunnel is set 


Toggle COARSE DOWN until water level is below the 10 cm 
mark then momentarily depress STOP (button with red light); | 
secure power switch; enter time into tunnel logbook 


EMERGENCY STOP: 


Secure switch on red box near breather slot 


B. CORRECTION FACTORS 

Data taken from within the wind tunnel must first be corrected and calibrated due 
to the physical geometry of this particular wind tunnel. Of primary importance are the 
actual Reynolds number and dynamic pressure being experienced by an object inside the 
test section. Since these are functions of velocity, an accurate method of obtaining test 


section velocity must be available. This velocity may be derived from the difference p, - 


Figure 5. N.P.S. Low Speed Wind Tunnel Water Manometer 

P,, but as was mentioned earlier, the tunnel Ap is meaningless unless it has first been 
calibrated to a known standard. This calibration procedure (which should be repeated 
before any important research is undertaken), is simply the tunnel calibration. The actual 
procedure and results are discussed in the next section. Here it is pointed out that 


Ap,once properly calibrated, may then be used to provide the test section freestream* 


*"Freestream" will henceforth refer to any parameter measured in the test section 
which is unaffected by any object (eg. an airfoil) placed in or any mechanism directed 
into (eg. a jet of air) the test section, and will be shown with (00) as the subscript. 


velocity V,., aS demonstrated by equation (1). | The development of the tunnel 


calibration equation is explained in the next section. 


A p = (calibration factor) x q = (calibration factor) V2 (1) 


Besides calibrating for q, there is a correction for the turbulence within the test 
section called the “Turbulence Factor," or 7F. It is found experimentally and was not 


repeated for this thesis [Ref. 1:pp. 32-33}. 


TF = (2) 


The TF was found experimentally for the NPS wind tunnel to be 1.04. The TF 
effectively raises the Reynolds number. 

The final corrections of interest to this lab are the tunnel boundary corrections due 
to solid and wake blockage caused by the presence of the cylinder in the test section. 
These corrections account for the test section walls and the apparent "closed system" the 
data acquisition system experiences. These are corrections for two-dimensional testing. 
Solid blockage correction €,, compensates for increased velocity and dynamic pressure 


due to the blockage caused by the model [Ref. 1:p. 36]. 


(3) 


K, is 0.52 for a model spanning the tunnel height and C,,, is the effective cross-sectional 
area of the test section. 

The second factor is é,,,, or the wake blockage correction. The decrease in velocity 
caused by a wake results in an increased velocity outside the wake as the "closed system" 
tries to maintain a constant mass flow rate within the test section. This results in an 


increased dynamic pressure on the model. [Ref. 1:pp. 36-37] 


(4) 


Here, d is the diameter of the cylinder, w is the width of the test section and C,, is the 
uncorrected drag coefficient. Reference 1 (pp. 36-37) goes into broader detail about 
what C,, represents. For the purposes of this experiment, C,, will be chosen to be c, 
(and thus labeled "c,,") found by each of two methods described in the next chapter. 
Reference 1 offers a simple estimate for total blockage on page 37, so no significant loss 


in precision will be experienced in calculating c,. 


Ca = Cg (1 - 3e,, - 2e,,) 


Thus (Ref. 1:p. 37], 


_ K,-Volume d 
i a A (6) 
2 
Cuz 


ia 1S then used in further corrections to velocity and the Reynolds number [Ref. 1:p. 


38]. 


V= V1 + ea) (7) 


Re,rpe = Re (1 + €,,..) (TF) (8) 


C. CALIBRATION AND RESULTS 
The relationship between Ap and tunnel q is essentially linear and provides the 
tunnel calibration curve. The slope and intercept of the calibration curve comprise the 
tunnel calibration factor which is the object of the calibration. 
The tunnel was calibrated on 26 April 1993 using three steps: 
1. The output from a digital manometer was compared to that of a known 
standard and a correction factor was determined for the digital manometer. 
2. The reading from the water micromanometer used for the wind tunnel (Ap) 
was compared to that of the digital manometer, which now read wind tunnel 
q using a portable pitot-static tube placed in the center of the test section. 


3. | The data obtained were reduced to a plot of Ap versus corrected tunnel q; 
this plot was the calibration curve sought. 


1] 


1. Calibration Step One 
_ The digital manometer was simply a transducer with one pressure port labeled 


“high” and one pressure port labeled “low.” The high port was connected to the output 
of a known standard, in this case a device called the "Schmidter." Next, the manometer 
reading of the Schmidter was compared to the output of the digital manometer over the 
pressure range necessary for the test. The "Schmidter" is shown in Figure 6, and 


includes a cylinder which can be manually pressurized and a water manometer reading 


in centimeters. 


Water (cm) 


Input a pressure and 
read 172 of scale 


Figure 6. The "Schmidter" (a calibration manometer) 


Two tests were made and the results were then averaged and plotted. Figure 7 shows 
the resulting conversion derived from this comparison. The subsequent slope and 


intercept (offset) of the plotted line were then applied to data found during step two. 
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Three things are to be noted conceming Figure 7: (1) the digital manometer was 
designed to read out in inches, not centimeters, and thus the conversion shown in 

Figure 7 shows inches; (2) it can be seen from the calibration equation in Figure 7 that 
the digital manometer had a linear error. This explains why the digital manometer had 
to first be calibrated to a known standard before it could subsequently be used to 
calibrate Ap; and (3) finally, since the digital manometer used had an operating limit near 
13 centimeters, higher runs were not pursued. Even though the calibration was not 
conducted at higher pressures, other data are available which ensure that the calibration 
factor remains linear for Ap’s above 13 centimeters. The AA2801 Cylinder Drag 


Experiment uses the wind tunnel in 4th gear at 15 cm H,O Ap. 


2. Calibration Step Two 

The next step involved placing a standard pitot-static tube into the wind 
tunnel, extending approximately 10.5 inches from the top center of the test section and 
facing into the freestream. The tunnel was operated in increments of | cm H,O (tunnel 
speed gauge (Ap)) and air velocity (in MPH) was recorded from the tunnel airspeed 
indicator. At each increment, the digital manometer output of the pitot-static 
tube was recorded (see Table 3). The tunnel was operated first in 4th gear, and then the 
operation was repeated in 3rd gear. The values read by the digital manometer were then 
averaged. The correction from Figure 7 was applied to this average and then plotted in 


Figure 8. 
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Correlation between the 

Digital Manometer and the 
Digital "Schmidter" (Howard/Miller 
Manometer 26 April 1993) 


(inches) 
5 


Yt ‘ 
M Actual infcm conversion 
“ye 


0123 45 67 8 9 10 11 12 
"Schmidter" (centimeters of water) 


Figure 7. Data From Step One of the Calibration 
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TABLE 3. Step Two Results 


hcpeusiNecadd Gani. | Icoreeeted q (cm) | 


3. Calibration Step Three 
Step three was to plot the first and fourth columns of data from Table 4 
against each other. Again, a simple curve fitting program was used to find the 
calibration curve. Here, the slope and offset of Figure 8 provide the tunnel calibration 
factor, and the equation is 0.243 + 0.895 * q, with q given in cm. of water. This 


differs from an earlier calibration of 0.93 * q. It should be noted that the two formulas 


15 


are merely curve fits to measured data, and do not vary as much as might be supposed. 
The curves produce identical values at about 7 cm., and at 15 cm. only differ by about 
0.25 cm. At 10 cm. the difference is about 0.1 cm., well within the fluctuation level of 


the micromanometer reading during the periodic pressure surges in the tunnel. 


NPS Tunnel Calibration Factor for the NPS Aerolab 
Low oredr Wind Tunnel Howard/Miller - 
(26/93...plot revised 12/3/93) 


Tunnel "Delta 
P* (cm H20) 


13 
12 


Tunnel “q" - ss measured by digital manometer 
and corrected usin 
q(cor)=(q(uncor)-0.04668)/0.41018 


Figure 8. Calibration Curve from Step Two data 


Dynamic pressure for incompressible flow is 


tunnel q = q. = 52 (9) 
and Ap (in lb/ft?) is 
Ap = [py (em) ~ Py (MI aud = ie (10) 
where 
: lb 
Ywoter ~ 02.35 eo (11) 


The wind tunnel calibration equation, from Figure 8, is 


Ap = 0.243 (cm H,O) + 0.895 9,05 section (12) 


And combining equations (9) through (12) and taking density p,. in slugs/ft’, 


Vv. (fy 2 7s - 0.243) 
sec 0.895 p. 


This equation has been written into a MATLAB code which uses different unit 


conversion factors to create updated velocity tables. This code is contained in Appendix 


E. 
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I. CYLINDER DRAG THEORY 


A. AERODYNAMIC DRAG CLASSIFICATION 


From Reference 2: 


TABLE 4. Drag Classification. 


ee 
| TOTAL DRAG - D | 
| 
| 


| SHEAR - D, NORMAL 
(due to pressure) 


(due to skin friction) 
; Form - D, Induced 
(due to separation) (due to lift) 


Dominates streamlined 
bodies 


Dominates blunt bodies 


Profile Drag = D, + D, (2-D flow) 


Parasite Drag = D, + D, (3-D flow) 


Though not used in this experiment, a wind tunnel balance may be used to directly 
measure the drag (and lift/side) forces on a cylinder in the test section. The balances 
used are divided into two categories: internal (sting-type) and external. The NPS Low 
Speed Wind Tunnel has provisions for an external balance using a strain-gage bridge 
network which measures lift and drag. 

Methods which use pressure to find drag include the pressure distribution method, 


which finds form drag only, and wake analysis, which finds profile drag. 
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B. MEASURING DRAG BY PRESSURE DISTRIBUTION 


1. Theory. 


0 al 


Radius =r 


Figure 9. Pressure Forces on the Surface of a Cylinder 


As shown in Figure 9, which was derived from a similar figure in Reference 3 (p. 
58), there is a pressure force df acting normal to the surface of an airfoil. There is 
another pressure force not shown which acts tangentially to the surface. This tangential 
force is called shear and is a result of skin friction. However, for blunt bodies, drag 
related to surface pressure dominates and the tangential pressure force is commonly 


ignored (Ref. 3:p. 60]. So, the force shown is the significant contributor and has a 
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component called form drag dd associated with the streamwise drag acting on the 
cylinder. This form drag may be found by starting with the normal pressure force as 


follows: 


df = p da = p (r d8) (14) 


where df is acting Over a unit span and r is the radius of the cylinder, as shown. It 


follows that 


dd = dfcos® = p-r-d8cos6 (15) 


where dd is the form drag component and di, shown in Figure 9, is the induced drag 
component. By integrating over the whole surface of the cylinder to sum all the dd 


components, 


2x 2n 
= 1 
: dd rf p cos6d6 (16) 


which equals the form drag, or 


Di=r f, *D cos8d8 (17) 
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The static pressure p acting normal to the cylinder’s surface can be 
considered as a local static pressure which varies with the pressure port location (as the 
cylinder is rotated from 6 = -20° to 6 = +180°) and the freestream static pressure 
which is constant. The freestream static pressure p,. can be pictured as a “blanket” of 
constant thickness enveloping all 360° of the cylinder. Integration of this “blanket" over 
a closed surface will result in zero. 

As noted previously, integration of the local static pressure over the cylinder acting 
in the downstream direction will produce the form drag. Calling this local varying static 


pressure Pryis 


2n 
D,=9r f, P,,£080d8 (18) 


And due to symmetry, 
[;"cosed8 = 0 (19) 
which can be written due to constant p., 
r [;"p.coseae = 0 (20) 


Therefore, the part of the pressure which contributes nothing to the form drag can 


be removed from the previous expression: 


2n 2n 
D, = rf Poy COS8d8 - rf, p.. cos6d8 (21) 
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which simplifies to 


= rf (Poy ~ Pe) cos8dd (22) 


Furthermore, since the drag equation is already per unit span, it can now be made 
dimensionless. Also, as two-dimensional flow is being examined, the total reference area 
S reduces to the chord length c x 1, which for the cylinder is 2r x 1 and 


2a 
Dy _ "fy Can ~ P=) costco 3) 


Cp = | Ga 


» q. § q.. 27 


The r’s cancel, leaving 


1 2 
a P.1 P.) cos6d6 (24) 


a = 


The pressure coefficient, C,, is a dimensionless ratio between the local pressure acting 
normal to a given surface minus the freestream static pressure (p,.), and the freestream 
dynamic pressure (q..). With reference to equation (24), since q,, at a fixed tunnel speed 


is simply a constant, 


—— Pe ) cos0d® (25) 


4, 


22 


which is the same as 


1 f2« 
4,2 2 
4, = af, Cp c0s0a8 (26) 


The cosine function is an even function and by symmetry the last equation can be written 


as 


The following authors arrive at the same conclusions using slightly differing 


approaches: 
1. Anderson, [Ref. 2:p. 209, equation 3.130) 
2. Bertin and Smith, [Ref. 3:p. 59, equations 2.40, 2.42] 
3. Rae and Pope, [Ref. 4: p.221, equation 4.29] 
4. Hoerner, [Ref. 5:p. 1-9, equations 4-6] 


There may be some concern regarding the sign convention used in this 
experiment as presented here. In the AA2801 Cylinder Drag Lab, it is assumed that 0 
= 0 and 6 = =x are located as shown in Figure 9. Thus, integrating an otherwise 
dominating positive area yields a correspondingly positive c,. Classic discussion 
normally has the limits of integration reversed, with 6 = 0 beginning at the leeward side 
of the cylinder and growing positive in a counterclockwise fashion. Additionally, after 
seeing a plot such as Figure 10, which compares theoretical (inviscid) C, to actual 


(viscous) C,, one might be led to believe that a negative sign would be required in front 


of Equation (27) just to make c,, positive. 


fron. fre abs amg eta 


The ree tac ald 


from Cp } 4sin ?tthetay [Andersen 


Figure 10. Inviscid C, as Compared to Actual C, 
Actially. the integrand is not simply C,. this C, ~ cost. Figure PP makes this apparent 


and addresses the positive -area/positive sign question. 


Cp only 


Figure U1. C,“cas# compared to C,, 


C. MEASURING DRAG BY WAKE ANALYSIS (MOMENTUM) 


1. Theory 
In Reference 6 (p. 66). steady 1-D flow (as ina wind tunnel test section with 


no flow distortion) can be represented by the momentum equation: 


Fe Cee a: (28) 


& 


c 


Drag. by definition, is a negative component of YF. so swapping V., and V,,. in 


equation (28) will keep their difference a positive value. The flow behind the cylinder, 


for purposes in this laboratory, is two-dimensional. Rae & Pope develop the momentum 
equation (Ref. 4:p. 214] into a double integral and state that “...the part of the air that 
passes over the model suffers a loss of momentum, and this loss is shown by and equal 
to the profile drag..." Profile drag, as mentioned earlier, is two-dimensional drag 
composed of both skin friction (D,) and form drag (D,), and is collectively denoted here 
as D. Thus, since mass flow rate equals pAV, where A is the area of the wake 


perpendicular to the freestream, 


D = [fevaav. -V) (29) 


In incompressible (low speed) theory, p, equals p,. (a constant) and the equation 


becomes 


D = [f(p.v.v, - pV) dA (30) 


For a unit span of the cylinder, as was explained earlier, S equals c x | equals 2r x 1, 
and for a unit slice of the wake along the y-axis, Da equals dy x 1. So two-dimensional 


drag, denoted by d may be written 


d= fIp.vV, - pv) dy x 1 (31) 


Since 


D 
Cc, = — 32 
D 7s (32) 
then 
4 afc x1) ( 1 oV.2)2r (p V2r (33) 
2 eo ae 
and substituting equations (31) into (33) yields 
V.-pV2 
~ - LbeVV “0M dy ai 


p.Vor 


The densities cancel out and since the rest of the denominator is also a constant, then 


VV. v2 
== ft=2-4y (35) 


Equation (9) from page 17 may be rewritten as 


v= | | white v, = | 2 (36) 
P. Pe 


Substituting equation (36) into equation (35) yields, 


ars ft |= - 2% an 


7 q. 


Since q,, and q.. are simply ((Op),, - Py) and ((Po). - P..), the difference between these two 


ratios are then integrated over the width of the wake to get 


Similar results to those in the previous analysis were obtained by: 


Anderson, (Ref. 2:p. 107, equation 2.74] 
Rae and Pope, [Ref. 4:p. 215, equation 4.24] 
Hoerner, [Ref. 5:p. 1-7] 

Panton, [Ref. 9:p. 386] 


PON 
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IV. AA2801 CYLINDER LAB 


A. PREVIOUS LABORATORY 


1. Procedures 
The former laboratory procedures for the AA2801 "Drag on a Cylinder" may 


be found in Reference 7 and are reproduced in Appendix A for contrast. 
2. Limitations and Problems 


a. Limited Data Acquisition 

Chief among the limitations of the previous lab was the lack of 
numerical data. Thus, all analysis was strictly drawn from graphs made by the Hewlett- 
Packard X-Y Plotters (see Figure 12). Students would count squares under the curves, 
assign coordinates to the curve, use a planimeter to measure the area under the curve, 
and one student was even known to cut out the curves and a small representative square 
and then weigh them on a precision digital scale! Also, not being able to retrieve data 
meant that this lab served only one real purpose: Estimating Drag on a Cylinder. 
Certainly the potential to expose students to the functional use of the wind tunnel existed 
but was not exploited by this approach. Beyond the lack of data retrievability, there 


were other problems with the former lab procedures. 
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Figure 12. Hewlett & Packard X-Y Plotters. 


b. Poor Relation Between Method Chosen and Actual Procedures 


The wake analysis equation given in the old lab, 


2 V, V,., 1 P, - P, 
C, = ~=f{— - (—) by - —/[——]6 (39) 
ea Gea 


was supposed to have been derived from Newton's Second Law using the Reynold’s [sic] 
Transport Theorem. To further complicate things, the lab had students allegedly plotting 
two (but not all three) parts of this equation. These two parts were: 
(V,/V,)? vs. y and (P, - P,)/q; vs. y. 
Next, the students were to subtract one curve from the other to estimate 
C, (not written as "c,"). To begin with, the equation does not agree with contemporary 
definitions of the wake analysis method (see Chapter III, Section C). Secondly, the 


procedure does not support the equation - either intuitively or precisely. While it was 
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pointed out that (V,/V,)? was simply q,/q,, the V,/V, portion was missing. Also, 
subtracting one curve from another does not lead to the solution of the equation given. 
Finally, the P, - P, values do not represent anything useful. P, was the static source in 
the wake while P, was the static port on the cylinder directed into the freestream at 0 
equals 0°. The difference, measured this way, was between the static pressure through 
the wake and the stagnation pressure at the cylinder. | This was confusing and 


bothersome for an otherwise simple method of drag estimation. 


c. Complexity of Procedures and Sparsity of Educational Value 

To make the plot described above, the student had to crank the traverse 
over and back while the technician had to change leads into the X-Y plotter. Graph 
paper had to be placed properly and the technician was inordinately relied upon to make 
the lab work. The student played a minor role in the lab other than observation and 
cranking. Perhaps one of the more valuable parts of the experiment, the setting of the 
scale of the transducers, was left to the technician. Additionally, the student leamed very 
little about data acquisition, something that would greatly benefit him or her for future 
use of the wind tunnel. The lab handout provided sparse information for answering the 


question, "Why?" 


d. _Inconsistency/Inaccuracy of Terminology 
P, was never defined, nor was q,. "C," was used for two-dimensional 
flow where "c," is the convention used by most texts. The length of the cylinder was 


given as 28.0 inches when it is actually 28.4 inches, which changes the volume from 
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192.68 in’ to 195.43 in’. The pitot probe and the static probe on the traverse were 


referred to.as "two pitot-static probes.“ 


B. UPGRADED LABORATORY 
Data Acquisition and System Description 


a. Schematic 


Figure 13 illustrates the signal path for the AA2801 Cylinder Drag Lab 


upgrade. Meanwhile, the photograph in Figure 14 is significant: data acquisition, data 


storage, data retrieval and data reduction - all in one place. 


Figure 13. Signal Path of Data Acquisition System 
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Figure 14. 486DX/33 MHz PC Used For Data Acquisition 
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b. Signal Configuration 


(1) Py - Po. There are two static pressure ports: a local pressure 
port (p,,,) located on the surface of the cylinder approximately mid-height and a static 
pressure port (p,,) located on a traverse probe mount. Both pressures feed into a Statham 
pressure transducer rated at +1 psid, 15 Vdc maximum input. The transducer is excited 
by approximately 11, volts and the output signal is then fed to an NPS-built signal 


conditioner, and then on to the PC. 


(2) (Dow ~ Py equals q,. Traverse mounted pitot and static probes 
((Po)w and p,,) are co-located behind the cylinder and the pressure measured from these 
are fed into a Statham transducer, Model P69-1D-350 rated at +1 psi. The signal path 


to transducer and to the PC is identical to that of the previous transducer. 


(3) y. # mechanical traverse operated by a hand-crank moves the 
mount which holds the pitot and static probes behind the cylinder. A small gear drives 
a 10-turn, SOK Q IRC potentiometer, Type H115-T-5620. The potentiometer is excited 
with approximately 11 volts and the output signal is fed from the NPS-built signal 
conditioner to the PC. 

(4) 06. The cylinder is mounted to a turntable, with rotating gear 
driving an Amphenol 100 Q, 10-tum potentiometer, Type 4301 B. The potentiometer 
is excited by 0.56 volts while the output signal is fed to a Pacific Amplifier/Signal 
Conditioner, Model 8255 and amplified 100 times before routing to PC. The location of 


the turntable’s azimuth scale is such that the cylinder local pressure port is pointing 
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directly into the freestream when 6 equals 0°. The turntable rotates the cylinder 


clockwise from @ = -20° to 6 = +180°. Figure 15 shows the rotatable base. 


Figure 15. Electrically Rotatable Base for Cylinder 


(5) (Pda - Po equals q,,. The tunnel airspeed indicator pitot-static 
tube, located under the ceiling in the forward part of the test section, is used to provide 
(Po). and p,, to a Statham pressure transducer rated at +1 psid, 15 Vdc maximum input. 
This transducer is excited by approximately 11 volts and returns its output to an NPS- 
built signal conditioner (shown in Figure 16), then to the PC. The airspeed indicator 
serves as a reference static pressure under the assumption that p,. read there remains 


essentially constant at a fixed Ap, as was verified by calibration. 
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Figure 16. N.P.S. Signal Conditioner 


(6) | CIO 37-pin MINITERM Connector Board. Al Signals arrive 
here first and are attached via "channel" points to a standard 37-pin RS-232 cable (see 


Figure 13). 


(7) ComputerBoards CIO-DAS16/330 A/D Card. <A high-speed 
analog and digital data acquisition card was placed in the PC with a base address of 
300H (300 HEX). A special library in the software package allowed the computer to use 


the A/D card (see Appendix B). 


(8) IBM 486DX/33 MHz PC. A standard 486 IBM-clone PC was the 
terminal for the data acquisition system. This PC had DOS version 5.0 and Windows 


version 3.1 loaded. 


c. Software 

LabTech Notebook for Windows, version 7.1.1, was the operating 
software used to drive the data acquisition system. It is a process and control type 
application which can be used to both collect data and to control systems. 

LabTech Notebook has two basic operating modes: BUILD-TIME and 
RUN-TIME. As their titles imply, BUILD-TIME is used to construct set-ups which are 
then executed in RUN-TIME. These set-ups contain customized "ICONblocks" which 
connect signals and create displays and data files. 

2. “ZERO" and "SPAN" 

These two LabTech Notebook set-ups were created to handle the job of 
scaling the signals coming into the program. It was felt that the signals must be correct, 
not just leaving the signal conditioners, but more importantly, entering the PC. Noise 
and other anomalies affect the signal at all points in its path. ZERO and SPAN are fully 
explained in Appendix B. 


3. Procedures 


a. Run #1 - Pressure Distribution Method 


(1) Approach. While maintaining the reference static pressure p, 
outside the wake and behind the cylinder (which is thus p,), the cylinder is rotated at a 
constant rate (and so is p,,,) from 0° to +180°. The turntable is started at 6 = -20° to 
ensure any acceleration is essentially zero by the time it passes 6 = 0° and the maximum 


P.y: iS actually recorded, although this is not required. 


(2) Measured sources. The signals being measured were p,,, - p, V. 
6, where p,, is the reference pressure and is assumed to be the freestream static pressure 
(p..) of the test section and where 6 is the position of p,,,, the local static pressure on the 
cylinder, as it varies from the stagnation pressure, through p,, past the point of 
maximum velocity, past the separation point of the boundary layer and finally to @ equals 
180°. 

(3) Purpose. These measurements will be used by the software to 


find and plot C, v. @. 


(4) Signal scaling criteria. Before running the tunnel, at zero tunnel 
Ap, since q, Pryi, Pe» Py and C, all equal 0, the transducer output is conditioned to also 
read 0.00 volts. Next, the tunnel is run at 15 cm. H,O, with @ equal to 0° and y equal 


to 0.0 inches. But at @ equals 0°, p,,, equals (p,)... Since it is known that 
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(Po) equals q.. + p., the following substitution can be made for C,: 


atd-= 0, <2 Se Fe Es 2 Be 4 (40) 


Knowing this, the output signal can be conditioned for the transducer to read 0.01 volts. 
The software used on the PC is then simply configured to multiply the signal coming in 
by 100 so that 0.01 volts is displayed as 1.0 when C, actually equals 1.0. Figure 17 


illustrates the LabTech Notebook set-up "PRE DIS!." Figure 18 shows the actual 


ry Nis 


3 
Sail 
3 
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Figure 17. LabTech Notebook Set-Up "PRE DIS!" 


results of executing PRE DIS! in RUN-TIME. 


(5S) Expected results. From inviscid, incompressible theory, the 
pressure distribution method will generate a cosine-type curve starting out with C, equal 
to 1.0 at @ equal to 0° and decreasing through C, equal to 0 to a minimum of C, equal 


to -3.0 or separation, whichever occurs first (see Figure 18). 


Figure 18. Actual LabTech Notebook Readout of C, vs. 6 from 
Run #1 


The point of separation will also define the point of wake initiation. It is anticipated that 
the wake behind a stationary cylinder in moving fluid, where the Reynolds numbers are 
of the order 2.7 x 10°, will not change appreciably in width (in the y-direction) from 
where it begins (Ref. 2:p. 238, Fig. 3.40d]. As the wake progresses in the x-direction 


away from the cylinder, it seeks to recover back to the freestream condition. Thus, the 
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width of the wake, for this Reynolds number, may be approximated by entering the value 


for azimuth at the point of separation into the following equation: 


maximum expected wake width (y..,,) = 2 (d | sin 8 | ) (41) 
where d is the diameter of the cylinder in inches. This value was used to predict the 


beginning of the wake during Run #2. 

(6) Using the data obtained during data acquisition to estimate the 
drag coefficient. Using the pressure distribution method for blunt bodies, form drag is 
the significant contributor to the resultant force on a non-rotating cylinder in airflow 
[Ref. 3:pp. 59, 60], [Ref. 2:p. 60, 209}: 


c, «= [* 775 cost 88, where C _ PP. (42) 
4% Jo q. > 4. 


The data obtain during Run #1 is a column of raw C, values, a column of corresponding 
azimuth (8) values, and a column of theoretical C, values, the latter of which are 
obtained from the equation for theoretical inviscid, incompressible surface pressure 


coefficients over a cylinder [Ref. 2:p. 199): 


C, = 1 - 4sin*e (43) 


Having these data, any program or technique of integrating the 
first two columns over the interval of 0 to x will supply the value c,,. Tabular data were 


not available in previous labs and students relied on physically measuring the area under 
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the curve from 0° to +180°. Now, using LabTech Notebook and automatic data 


acquisition, data are stored in an ASCII file and may be imported into a program such 
as MATLAB for analysis. 
To compare the c,, obtained from the experiment to empirical 


data, the Reynolds number is calculated from 


Ree (44) 


and the graph of Re vs. C, shown in Reference 7, (p. 36) or Figure 24 is used. Care 


must be taken to correct Re for the turbulence factor, TF, to produce the effective Re: 


Reger = Rsncorrected.t +) ‘TF (45) 


where 


et, te. (46) 


b. Run #2 - Wake Analysis Method 


(1) Approach. While the cylinder is fixed in some (any) position 0, 


the wake probes are moved along a y-axis traverse from 0.0 inches to +20.0 inches. 


(2) Measured sources. Three values are measured during this run: 
Gu» Ge, and traverse location, y. By definition, (Po),, - P. equals q,, and (po). - Po equals 
qe. It has been shown experimentally that when the probes which measure (Po), - Py are 
located outside the wake, q, equals q.. Since q,/q.. is the ratio sought, the wake 
analysis (momentum) method is really used to determine variations in velocity through 


the wake. 


(3) Purpose. These values will be used by the software to plot the 


following difference in ratios: 


Pow = Pw _ Pow = Pw yy (47) 
(Pe - Pu (Pi). - Pe 
which is 
fe _ ae yy (48) 


The difference between these two ratios forms the integrand of the momentum method 


for calculating the drag coefficient, as was explained earlier. 


(4) Signal scaling criteria. Before running the tunnel, at zero tunnel 
delta P, all pressures are equal, all q’s are zero and all transducer outputs are scaled to 


read 0.0 volts. Next, the tunnel is run at 15 cm H,O, with y equal to 0 inches. Here, 


outside the wake, q,, equals q,, (both # 0.0). Since q,/q.. will be plotted, using the 


Commutative Law, both outputs are scaled to the same value, in this case 1.0. Actually, 
any common value would have worked since the ratio when q, equals q.. is always unity. 


Figure 19 shows the set-up WAKE ANI/2 from LabTech Notebook. 


Hake tnd 2 Set-up 
Lablech Netebacdk 
Mutdows ms 
Vey id | { 


Figure 19. LabTech Notebook Set-Up "WAKE AN1I/2" 


(5) Expected results. From the value of maximum expected wake 
width (yy) estimated during Run #1 (equation (43)), no appreciable change in pressure 


should be seen until reaching the interval 


/ 
Yware (inches) = 10 + wa (49) 


Logically, q,. should remain constant during the tunnel run. However, q,, the local 
dynamic pressure in the wake, should change to reflect the condition of the wake as the 
traverse moves from 0.0 to 20 inches. Also, p, should change less than (p,),. Indeed, 
many "profile drag” rakes used to measure wake drag from attached flow over airfoils, 
often have only a single static probe for dozens of total pressure probes [Ref. 4:p. 112). 
The static pressure is treated to essentially remain constant through the wake. Actually, 
due to tunnel wall effects, the static pressure may change slightly [Ref. 1:p. 35). 
Nevertheless, it is expected that as the wake is entered, q, should get smaller because 


(Po), decreases at a faster rate than p,. 


(6) Using the data obtained during data acquisition to estimate the 
drag coefficient. The software set-up "WAKE ANI" in LabTech Notebook is 
configured to receive the three inputs: q,, q.. and y. Through a string of "calculation 


blocks," WAKE AN] forms the plot mentioned earlier of 


Pow ~ Pw Ody ~ Pw 
(Po). - Pu (Py). ~ Po 


vy 


which is 


The ordinate of the plot is the integrand from 


1 5%, | PodwPw _ Podw-Pw (52) 
Ca = ( Oi. ra a 


This equation is integrated from 0.0 to 20.0 inches and the result is an approximation of 
the drag coefficient. Although the setup WAKE_ANI creates a display of the ratios of 
dynamic pressure versus traverse distance y, three columns of data are also generated for 
the student to analyze, using a program such as MATLAB. They are simply the three 
inputs q,, q. and y. It is a simple matter to write a code to use these inputs to solve 


for cy. See Appendix D for an example. 


c. Run #3 - Wake Analysis Method with Tripped Boundary-Layer 

The only difference between Run #2 and Run #3 is the addition of two 
"cellophane" adhesive tape strips, 4" x 20" x 1 mil. These strips are affixed to the 
cylinder, approximately 1.25" on either side of the @ equals 0° position, parallel to the 
z-axis (vertical axis). 

The function of these two strips are to act like tiny “vortex generators," 
which are familiar to all Naval aviators. The strips effectively "trip" the attached 
boundary layer causing otherwise laminar flow to become turbulent. The consequence 


of this act is dramatic. By inducing an early transition from laminar to turbulent flow, 


the engineer can delay actual flow separation, thus delay the formation of the wake. This 


in turn decreases the drag produced by an object in an airflow, a favorable situation for 


blunt bodies [Ref. 2:p. 644]. 

There is an intuitive explanation for why flow separation is delayed by 
creating turbulence where there once ea laminar flow: In the very small boundary 
layer, laminar flow can be thought of as flow in 1 dimension only - parallel to the 
surface. There are no significant forces at work to move the flow either closer toward 
the surface or farther away. Eventually, however, the displacement of the flow from its 
original freestream position becomes so great that the laminar boundary layer flow leaves 
the surface of the airfoil and seeks to return to its former place. If, however, this 
laminar flow is artificially brought into contact with the surface of the airfoil, which is 
by no means smooth, the air molecules will be caused to trip and tumble. Specifically, 
in some random direction in relation to the surface of the airfoil, but in general in a 
continued downstream direction. As this air trips and tumbles (turbulence) and moves 
along the surface of the airfoil, it acts like a cartoon snowball rolling down a winter 
hillside. As it grows, it "reaches" farther above the surface of the airfoil and “gathers” 
in higher energy air. This higher energy causes the air near the surface to become much 
faster than it was when it was laminar. This new speed in turn causes the flow to “push” 
ahead farther along the surface than it naturally wants to go. Eventually the forces at 
work trying to restore all flows to the original freestream condition overcome the contact 
between the now turbulent flow and the surface and the airflow separates and is returned 
to the freestream. Now the point of separation is much farther along and thus the wake 


width and associated drag is significantly reduced. 
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Anderson, in Reference 2 (p. 643), states: 


Because of the agitated motion in a turbulent flow, the higher-energy fluid 
elements from the outer regions of the flow are pumped close to the surface. 
Hence, the average flow velocity near a solid surface is larger for a turbulent 
flow in comparison with laminar flow...because the energy of the fluid 
elements close to the surface is larger in a turbulent flow, a turbulent flow 
does not separate from the surface as readily as a laminar flow. If flow over 
a body is turbulent, it is less likely to separate from the body surface, and if 
flow separation does occur, the separated region will be smaller. As a result, 
the pressure drag due to flow separation, D, [form drag], will be smaller... 


Though not demonstrated in this lab, Hoemer points out in his classic 
Fluid Dynamic Drag (p. 3-27) yet another way to reduce separation around a stationary 
cylinder: use suction on the cylinder via a porous surface. This prevents or delays the 
start of an alternating vortex sheet. 

The helicopter industry is already producing aircraft which use a 
delayed-separation principle called the "Coanda"-effect. The application is NOTAR’, 
which stands for "No-Tail-Rotor" [Ref. 10:p. 1f]. Here blowing, vice suction, is used 
to manipulate the point of separation around a cylindrical tail boom to provide an 


asymmetric surface pressure distribution, which then generates desired lift. 
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V. RESULTS 


A. ACTUAL CLASS DATA (USING DRAG1.M) 


Table 5 is a summary of the results obtained by using the program drag! .m created 
in Matlab code, which is located and explained in Appendix D. Users of drag].m may 
access the explanation by typing 

> >help dragl 

Dragl.m assumes the following: NPS Low Speed Wind Tunnel and associated TF 
and calibration factor; tunnel operated at 15 cm H,O; equations contained in this thesis 
are used; data have been edited to remove "idle data." 

Table 5 results are from actual classroom data obtained from the three LabTech 
Notebook set-ups described in Chapter IV. Two AA2801 labs of approximately seven 
students each were conducted using the upgraded laboratory handout contained in 
Appendix B. Since this lab was conducted on 9 November 1993, Appendix B has been 
revised to include the lab times required, Table 1 - Drag Classification, and illustrations 


under the ZERO/SPAN paragraph. 
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TABLE 5. DRAGI.m results 


» dragl 

% (from data taken during 9 Nov 1993 % 
AA2801 Lab) 

% 

% This data has been reduced using 

% DRAGI1.m. It has been corrected 

% for blockage and is explained in 


% 

0.71 

0.68 

0.67 

7.3 in 
12.7 in 
64.1° 

0.68 
2.89e+05 
0.99 


—_ 
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= 
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B. PRESSURE DISTRIBUTION METHOD (RUN #1) 


1. Data 

Data were obtained, as explained in Appendix B, and filed under 
PRE DISI.m. Actual results from the classroom are recorded in Appendix C. It should 
be noted that, in an effort to save space, this DOS Text File has been saved in small 
print and placed into multiple columns. 

The set-up PRE_DIS1 automatically begins collecting data when the rotator 
passes -17° and stops when it passes +179°. This worked because the rotator itself 
turned at a constant rate and was controllable electrically. The data were originally 


saved in two columns: @ and C,(raw). All data prior to 0° were manually edited out. 
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Due to the assumed constant rate of the rotator combined with the constant 


sampling rate of 10 Hz, the trapezoidal rule was applied to the data using [Ref. 8:p. 


193): 


ve pace > - ypl xf = j (53) 


An uncorrected (for blockage factors) c,, of 0.71 was obtained using a 


modified trapezoidal rule: 


. Ya-1) Ne 
A = Dak 2 2) x IX —Xq-ry! (54) 


Here, equispaced x-values were not assumed. 

The blockage factors described earlier (equations (3)-(5)) were applied to the 
results from both methods, yielding a c,, of 0.68 using equation (55) and 0.67 using 
equation (56). Figure 20 shows how the MATLAB program Drag! .m plots C, vs. @, and 
places an asterisk where separation is believed to occur. A routine which assumed C, 
increases abruptly during classical cylinder flow separation predicted the separation 
azimuth to be 64.1°. This value for theta was then used with equations (43) and (51) to 
estimate the wake width, which ranged from 7.3 to 12.7 inches. 

Next, Matlab-defined functions "polyfit" and "polyval" were used to generate 
a fifth-order polynomial approximation to the actual data. In this case, a fifth-order 


polynomial was chosen as it best fit the plot given. Results are shown in Figure 21. 
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Late from relish m 
KUN #1 


Pressure Ustrrbution 


Figure 20. The Separation Point Estimated by SEPAZ 


Finally, the effective Reynolds number of 2.89 x 10° was calculated using 
equations (7), (8), (13) and (46). 
2: Analysis 
The corrected c,,'s plot close to but not on the empirical data shown in 
Figure 24. One reason why 0.67 is seemingly low is that the pressure distribution 
method does not measure cy. it only measures the component due to pressure (form 


drag). In fact. Panton, in Incompressible Flow [Ref. 9:p. 386]. was responsible 
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Cp*cos(theta) vs theta 


Lels from Pre.Ds/m 
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Sth order potynonial 


C. WAKE ANALYSIS METHOD (RUNS 42,3) 


1. Data 
Essentially the same methods just explained were used to find c, related to 
the wake analysis from Chapter V, Section B. Equation (56) was used to find c, for 
Run #2 (0.99) and Run #3 (0.18). The wake analysis method uses the WAKE ANI/2 
set-ups and from data reduction using DRAGI.m, the following plots were generated 


(Figures 22 and 23): 


Nila from Boke Aplin 
KUN BE 
Hike Apofpays 


Figure 22. WAKE_ANI Results (Run #2) 
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Figure 23) WAKE AN? Results Compared (Run 73) 


Poouse these data (afse recorded in Appendix ©). the dbo pomts’ were tie 
edited ont Untike the fist setup during which the sampling rate was a steady LOH, 
the movement of the pitot static probes an the traverse was cntirefy dependent on the 
speed and accelenions of the studeat dotng the cainking  AddtGenatly . no pros istan 
wascmade fo automatically starbor stop the data acquisition A keypress imitated the dite 


collection and selecting FXNEP stopped it Consequently. idle data existed oo both ends 


of the data files. These idle data were easy to detect. All data prior to the y-value (@) 
increases and all data after it stops increasing were edited out. 


Figure 24 shows the results of all three runs plotted against emperical data. 


FIGURE 3.39 


Variation of cylinder-drag cocfiicient with Reynoids sumber. ( Experimental date as compiled in Panven, 
Ronald, incompressible Flow, Wiley- Interscience, New Yerk, 1964.) 


Figure 24. Results Plotted Against Emperical Data 


2. Analysis 

Again, the corrected c, from Run #2 was plotted against the Re,, and 
compared to the empirical data. It compared very favorably, tending to validate the TF 
given. Also, from Figure 22, the wake analysis width prediction of DRAG1.m seemed 
reasonable. 

The 81.7% reduction in c, from Run #3 is dramatic. If this number is high, 
then the TF might come into question again. "The higher the freestream turbulence, the 
more readily transition takes place [Ref 2:p. 389]." Golf balls use dimples to accomplish 


similar levels of dramatic drag reduction. 


There is one concern that arises out of Figures 22 and 23. It would appear 
that given the use of equation (53) it should be quite impossible to ever have a negative 
solution, as was seen during both runs. There are two possible explanations for this 
phenomenon: (1) As shown in Figure 25 and as discussed on page 35 of Reference |, it 
is theoretically possible to have q,, become larger than q,.. Were this to happen, the plot 
of equation (53) would indeed be negative for that interval. (2) A more likely reason, 
though, is that the signal from the q,, or q,. transducer is either poorly scaled or actually 
drifts with either time or with displacement. This could be easily verified by repeating 
a traverse of the test section in both directions, with the cylinder removed, and 
comparing q, to q.. Unfortunately, time did not allow for this verification and a 
solution to the negative values seen in Figures 22 and 23 are left for future exploration. 
For this thesis, the data were neither adjusted up nor down (from the zero axis), nor was 
the plot in Figure 22 pivoted by anchoring the left end and rotating the right end of the 
plot up to the zero axis. These corrective methods might result in a better 
approximation, but it remains to be seen whether data acquisition or tunnel velocity 


profiles are responsible for the negative values. 
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+ 


(see Ref. 9:p. 388) 


This plot shows then, =. (see Ref. trp. 35) 
what would happen if . 
q(wake) => q(infinity), that ts, 
if the ratio (qw’ql) was ever 
> 1.0 (shown here vielding a 
negative result on the interval 


Figure 258. Graph of the Equation y =sqrt(x) x 


VI. CONCLUSIONS 


A. UPGRADE LIMIPPFATIONS 


{. Wake Analysis Method 
Rae and Pope point out in Reference 4 (pp. 214-217), that the wake analysis 
method has several "pray areas:" 
I. How do you measure spanwise variation of the profile drag 


behind an airfoil? [p. 213] 
2. Are vortices present at higher Re? [p. 214] 
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3. Is true two-dimensional flow a function of Re or does it really 
exist? [p. 214] 


The authors state twice: 
It is well known that the wake survey is not valid...where separation is 
present...cannot measure momentum loss caused by separation or fluid 
rotation. [p. 214,217} 
The limitations, then, are that it is yet unknown exactly how valid the 
approximation of c, is from the wake analysis method. Certainly flow separation is 


present. While the value obtained, c, = 0.99, fits that empirical data, enough questions 


are raised by Rae and Pope that no strong conclusions may be drawn at this time. 


2. Hardware 
Current hardware returns signals on the order of 10° volts. Subsequent 
amplification of these tiny signals has the effect of amplifying noise or simply limiting 
the usefulness of the signal. The Analog-to-Digital data acquisition card used in this lab 
has the ability to resolve signals down to 4095 parts. In some cases, the only signals 
possible were confined to a range of 0.000 to 0.009 volts (see Table 3, Appendix B). 
Meanwhile, the smallest range LabTech Notebook can program down to read is + 0.625 


volts. This means that 


10.625 - (_- 0.625) volts = + 0.00033 volts resolution (55) 
(4096 -~ 1 ) counts count 


Dividing an actual useful range of 0.009 volts by this resolution yields an approximate 


effective resolution down to just 28 signals discriminating between 0.000 and 0.009 volts. 
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Put another way, while the A/D card can resolve down to 4095 parts, only 28 of these 


parts are being used, less than 1% of the card’s capability. 


3. LabTech Notebook 
This software application requires considerable time to learn. Having a 
qualified user available helps greatly. The lengthy tutorial, once mastered, provides 
an excellent foundation for using this software package to do data acquisition. 
Beyond the time it took to learn, it has three other limitations: 
1. apparent limited sampling rate of 1000 Hz 


2. very limited graphics capability 
3. high memory usage in the Windows environment. 


B. RECOMMENDATIONS 
1. Study the Wake Analysis Method Further 

A survey of the static pressure in the wake behind the cylinder should be 
taken and compared to the freestream static pressure. If the static pressures vary 
significantly, the cylinder should be replaced and the tunnel velocity profile measured 
along the traverse, paying close attention to the ends. The velocity along the traverse 
should be checked for speeds which are greater than freestream. Rae and Pope wam 
that q must be taken as far behind the cylinder so that p, equals p,. [Ref. 4:p. 216]. If 
the wake velocity remains less than the freestream velocity, check for signal noise or 


other faulty modes of the data acquisition system. 


2. Match Hardware to Software Capability 
_The 12-bit resolution of the A/D card as just explained, combined with as 
large an output as possible, would yield a higher resolution than currently available. 
This means that the signals into the PC should have as wide a range as possible, not the 
current 9 millivolt range output currently seen. This output is both a function of the 
signal transducer and potentiometer and signal conditioner. If the hardware used to 
acquire signal sources must remain, then the signal conditioner must be reconfigured to 
amplify the output. If the sources can be measured using devices which discriminate 
well below .033 millivolts level, then the signal conditioners need not be changed. 
3. Investigate Better Data Acquisition Software 
Investigate non-Windows applications for better response to memory 
limitations of PC. Also, recommend obtaining, for long-time use, a data acquisition 
program which has a broader graphics package than LabTech Notebook currently 
employs. 
4. Validate TF 
Since no attempt was made to validate the turbulence factor, the TF should 
be recalibrated using procedures discussed in Reference 1. TF is important and 
contributes approximately equally with the blockage factors in the calculation of the 
effective Reynolds number (equation (8)). 
5. Upgrade the NPS Low Speed Wind Tunnel deck 
The next logical step in the wind tunnel upgrade is to match the data 
acquisition capability with the actual layout of the tunnel test facilities. Layers of old 
thesis projects, built upon each other over the years, can effectively be replaced with this 


data acquisition system. Former procedures and techniques may be incorporated into the 
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set-ups in LabTech Notebook on an as-needed basis. Furthermore, the current collection 


of signal conditioners, signal amplifiers, multiplexer and wiring inhibits easy access to 
the test section. The PC used in this thesis must be located a full twenty feet away from 
the test section making lab observation and display observation all but exclusive events. 
A permanent wood lab table extending from the test section along the diffuser, would be 
a practical way to provide this easy access and give a place for lab equipment, in 
component form, to be set up. Furthermore, a system of annually justifying the 
continued presence of former thesis projects is required to prevent the problem of 
layering. Copies of thesis reports should be readily available near each project left intact 


after the author, student or professor, transfers from the school. 
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APPENDIX A: PREVIOUS LAB HANDOUT 


Drag on a Cylinder 


Objectives 


The drag on a circular cylinder which spans the height of the wind tunnel test section 
will be determined using wake analysis and pressure distribution methods, demonstrating 
basic wind tunnel methodology, including operating procedures, Reynolds number 
calculation and application of wind tunnel correction factors. 


Apparatus 


The NPS low speed wind tunnel is described in detail in the Manual for 
Low Speed Wind Tunnel Testing. A circular cylinder (length = 28.0 in., diameter = 
2.96 in., volume = 192.68 in’) is placed in the wind tunnel such that it spans the test 
section height. A single static port is built into the mid-portion of the cylinder. The 
cylinder itself, is attached to an electrically operated rotatable base. In the aft portion 
of the test section are placed two pitot-static probes, aligned such that one provides total 
pressure and the other static pressure, in the same longitudinal plane. These probes are 
attached to a traverse mechanism which enables lateral positioning (across the width of 
the test section) to be varied. Pitot-static pressures are transduced by strain gauges and 
output voltages fed to x-y plotters (y-information). Additionally, sensors on both the 
rotatable base and traverse mechanism provide position signals to the plotters (x- 
information). Thus, wind tunnel output consists of five electrical signals: 


Pree = P 

Wwake = 2 

Poytinder = Po 
Cylinder position (@) 
Probe location 


The tunnel will be run at a pressure differential of 15 cm H,O (approximately 115 kts.) 
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Wake Analysis Methodology 


Applying the Reynold's [sic] Transport Theorem to Newton’s Second Law, 


OF = o (momentum) 
bt 


and writing the resultant equation for a control volume (see Zucker: Fundamentals of Gas 
Dynamics, p. 65-69) results in: 


2/,"2 V, 2 17,P. - Py 
Cc = = r, es = ees r,] (4.1) 
AG “Gye [tay 


Pressure Distribution Methodology 


The forces on a cylinder may be modeled as shown in Figure 21. 


Neglecting the shear 
force: 


Pretsurd — bieas dl = phud 1 hs phe sad 


Ue dd = -pAdé cos @ 


Figure 21. Aerodynamic Forces on a Cylinder 


Below the critical Reynold’s [sic] number, the form drag is the dominant component of 
the resultant force on the cylinder. (Typically, total parasite drag is composed of 
approximately 80% form drag and 20% skin friction drag.) For blunt bodies (see Bertin 
and Smith: Aerodynamics for Engineers, p. 88-91), form drag can be approximated by 
the expression: 


C,~[,"C,c0s880 


ah, oe P 
C= { —§_ 1695688 (4.2) 
a $y q 


X-Y Plotter Calibration 


X-Y plotters can be calibrated to provide the following ratios: 
q 
; = (VJV,P vs. y 


1 


P, - P, 
q 


vs. y 


P, - P, 
q; 


From the curves of these ratios, the integrals of equations (4.1) and (4.2) may be 
evaluated. Calibration of the plotters is accomplished as follows: 


° (VV, vs. y: 


The plotter pen is positioned at the zero point prior to starting the tunnel (V, = 0). 
After the tunnel is at full speed, with probes located away from the wake (V, = 
V,), the plotter pen is positioned at the (V,/V,)? = 1.0, y=0 point. 

3. | The probes are then traversed across the test section (20 in.) and appropriate ratios 
are plotted. 


NO — 
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After the tunnel is at full speed, with probes located away from the wake (P, = 
P,), the plotter pen is positioned at the C, = 0.0, y=0 point. 


2. The probes are traversed across the test section and the appropriate ratios versus 
y are plotted simultaneously. 
P, - P 
8 1 vs. 0: 
q 

1. The plotter pen is positioned at the C, = 0 point prior to starting the tunnel. 

2. | With the tunnel running, the cylinder is rotated such that 6 = 0°, which puts the 
static port at the stagnation point (P, = P, + q,, and the ratio equals unity). The 
plotter pen is positioned at the C, = 1.0, 6 = 0° point. 

3. The cylinder is rotated to the @ = -20° position and plotting begins (from @ = -20° 
to 180°). 

Laboratory Procedures 

1. Determine atmospheric pressure: 

Pim = (baro. reading)+(temp corr.)+ (lat. corr.) 
(Latitude correction = -.0245) 

2. Determine temperature. 

3. Make required wind tunnel runs. 

4. Trip the boundary layer and make an additional run in order to determine the 


reduction in wake size afforded by a turbulent boundary layer. (Wake analysis 
methodology only.) 


Conversions and Constants 


aig = 3.8 x 10? Ibf-s/ft? 
R,, = 53.3 ft-lbf/lbm-°R 
°F + 460 


1 
1 psi = 2.036 in Hg 
q = AP/.93 (Calibration factor) 


Report Requirements 


1. Calculate the effective Reynold’s [sic] number and determine the corresponding 
empirical drag coefficient from Figure 22. 


2. Utilizing x-y plotter curves, evaluate the integrals of equations (4.1) and (4.2) and 
determine experimental drag coefficients for each method, applying wake and solid 
blockage correction factors. 


3. Determine the percentage reduction in drag afforded by a turbulent boundary layer. 


| 


4. Briefly discuss results. 


° 
? 
vu 
w VU. 
g 


Figure 22. Cylinder Drag Coefficient 
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APPENDIX B: UPGRADED LAB HANDOUT 


CYLINDER DRAG EXPERIMENT 


Class time: :50 

Lab time: 1:30 

Report time: 3:00 

Total time: 5:20 

***Before the Lab*** 

You will need to read this handout as well as pages 11-15, 32-38 of the NPS Laboratory 
Win nnel Testing. You will also need to bring either 3/2" 

or 514" formated floppy disk with at least 100 KB space available. 


Objectives 


A current method of PC-based automated daia acquisition using a combination of 
pressure transducers, potentiometers, signal conditioners and amplifiers will be used to 
collect and analyze data from flow over a non-rotating cylinder. Two techniques of drag 
determination (pressure distribution and wake analysis) will be used, demonstrating basic 
wind-tunnel methodology. In addition to these two primary objectives, students will be 
introduced to the NPS low speed wind tunnel operating procedures to include tunnel 
calibration and tunnel correction factor. 


Apparatus 
The NPS low speed wind tunnel is described in detail in the Manual for 
Low Speed Wind Tunnel Testing. Figure 21 shows the signal path for the data 


acquisition system. The dimensions of the cylinder are d=2.96 inches, h=28.4 inches. 
The dimensions of the test section are w=45.0 inches, h(above reflection plate) =28.4 
inches. There is a single pressure port near the center of the cylinder. The cylinder is 
mounted to an electrically rotatable base which is able to rotate from -20 degrees to 
+180 degrees in a clockwise fashion. Behind the cylinder, on a manual traverse, is a 
combination pitot-static probe apparatus. The probes may be moved along the traverse 
in the y-direction from 0 to +20 inches. The transducers, potentiometers, signal 
conditioners and amplifier which connect these measuring devices to the data acquisition 
card inside the PC are shown in Figure 21. Specific information about each device is 
located in Table 1. The wind tunnel will be run in 4th gear at 15 cm H,O Ap. This 


68 


"Ap" is a 15 cm H,O difference between the local static pressure entering the contraction 
cone (p,) and the local static pressure leaving the contraction cone and entering the test 
section (p,). Previous wind tunnel calibrations have experimentally shown this particular 
difference of 15 cm H,O to equate to approximately 115 MPH. 


In general, the calibrated NPS low speed academic wind tunnel equation is 


Ap = 0.243 (cm H,O) + 0.895 4.0 secsion (1) 
Dynamic pressure is 
1 
jaction q = sh Mn (2) 


where 
cee ee tb) ;_\h 3 
Ap Pry [P, -P2](CM)LY ase Ge (5 a ] (3) 
and 
Your = 62:35 (4) 


Combining equations (1) through (4) and taking density in slugs/ft’, 


vy. (fy - | 2(Ap - 0.243) 
sec ‘ 0.895 p_ 
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Pacific Amp 


Rotator _ 


Software used: 


LabTech Notebook 


| 37-pin MINTERM 
i__| 


Figure 21. Signal Path of Data Acquisition System 


Terminology 

Po = total (stagnation) pressure 

Pp = static pressure 

q = dynamic pressure = py - p 

(Po) = freestream total pressure (associated with freestream velocity) 
Pp. = freestream static pressure (associated with freestream velocity) 
q. = freestream dynamic pressure = (Pp). - Po 

(Pow = wake pressure measured behind the cylinder at the probe 


((Po)~ assumed to = (py),.. when taken outside wake) 
Static pressure measured behind the cylinder at the probe 
(p,, assumed to = p,, when taken outside wake) 


no] 
2 
Il 
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qv. = wake dynamic pressure = (po), - Py 
P.y: = pressure measured at the cylinder 
D, = Form Drag (due to separated flow) 


D, = Shear Drag (due to skin friction) 


C, = pressure coefficient 

Ca, = uncorrected drag coefficient 

Cy = drag Coefficient relating to total drag or profile drag 
Ca = drag coefficient due to D, only 

r = radius of cylinder 

6 = azimuth of the rotator (signal from potentiometer) 

y = position of the traverse (signal from potentiometer) 
Conversion Constants 


latitude correction = -0.0245 in. Hg 

temperature correction = (see Table 4-2 on bulletin board near tunnel) 
Beir = 3.8 x 107 (Ibf-s)/ft? 

R,, = 53.3 (ft-Ibf)/(Ibm-°R) 

°R = °F + 460 

1 ft = 30.48 cm 

1 cm H,O = 0.0142 psi 

1 psi = 2.036 inches Hg 

heat section = (Ap - 0.243)/.895 (NPS wind tunnel calibration equation) 
TF = 1.04 (NPS wind tunnel turbulence factor) 


TABLE 1. Drag Classification 


TOTAL DRAG - D 


SHEAR - D, NORMAL 
(due to skin friction) (due to pressure) 


Dominates streamlined 
bodies 


Induced 
(due to lift) 


Hee =D.+D,GD) 
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Pressure Distribution Method 

For blunt objects in subsonic airspeeds, form drag dominates the resultant force on a 
non-rotating cylinder in airflow. The pressure distribution method only measures D,. 
Using Bertin and Smith’s Aerodynamics for Engineers, p. 60 or Anderson's 
Fundamentals of Aerodynamics, p. 209: 


Jf > —— cos6 80, where 


Wake Analysis Method 

Comparatively, the wake analysis method measures not force, but a change in velocity. 
It is a ratio of dynamic pressures, which reduce to a ratio of velocities according to 
equation (2). The wake analysis method accounts for profile drag, or D, + D,. From 


Rae & Pope’s Low-Speed Wind Tunnel Testing, p. 215: 


- 1 iy Ow Pw = Cdw-Pw ees 


(P)).-P. oP 


Data Acquisition 


The heart of any PC-based data acquisition system is the data acquisition A/D 
(Analog/Digital) board. This device takes analog signals such as voltages and converts 
them to digital signals for the computer software to use. In the case of this laboratory, 
pressure signals and position signals, in the form of output voltages (Vdc) are first 
conditioned for optimum signal gain, then channeled to the A/D board via a standard RS- 
232 terminal. The A/D board simply plugs into the inside of the PC and has two settings 
which are of interest to the user: 


1. 8/16 Channel Select - 8-channel differential input (normal or floating) or 16- 
channel single-ended input. This lab uses 8-channel floating differential input (Channels 
0-7). The term "floating" merely refers to the method of providing a common ground 
to eliminate unnecessary noise. 
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2. Gain/Range - manual or programmable. The A/D board used in this lab has a 
programmable gain/range accessible through the software. The gain/range may be set 
differently for each input device and for all uses of that input device. Bipolar and 
unipolar refer to the polarity of the range. An example of bipolar gain/range would be 
+5V. Unipolar might be 0-SV. 


Once the signals have been properly fed into the A/D board it’s up to the software 
to use the signals to provide actual data acquisition. This includes accepting the data, 
displaying the data and storing the data. Data that have been stored (for example a file 
of two columns of data: elapsed time and temperature) may then be retrieved into a 
spreadsheet and plotting program (such as Lotus 1-2-3°) or into a math processor (such 
as MATLAB’) to be analyzed. The convenience of the PC allows the operator to do this 
at one place and at essentially one time. 


A typical method of calling a data file into MATLAB’, for instance, would be: 


> >load A:\pre_diel.m 


> > whos % will display size of matrix [m,n] (you will need m) 

> >x=pre_disi(1:m,1); % creates s column vector x of the first column of data in pre dis! 

> >y=pre_disl(1:m,2); % creates s column vector y of the second column of data in pre_dis! 
> > plot(x,y).grid % creates a plot of x vs. y 


In this AA2801 lab, the following data will be generated by LabTech Notebook: 


TABLE 2. File Organization 


po Columns | Column 2 | Column 3 


| PreDistim | @ | Cp Raw) | Cp (theoretical 
_ _ 
_ _ 
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If the A/D board is the “heart” of data acquisition then the software is the “brains.” 
Data acquisition software programs (such as LabTech Notebook’) allow the operator to 
configure “set-ups” which tell the PC how to collect, display and store the data. Current 
software uses block icons which may be cascaded into an infinite number of combinations 
of configurations. 


Safety Procedures 


Prior to executing the lab, normal pre-operating safety checks should be made. Students 
should use proper hearing protection for tunnel speeds above 5 cm H,O. 


Initial Procedures Prior to Running the Tunnel 


1. 


Determine atmospheric pressure (p,,, = (baro)+(temp. corr.)+(lat. corr.)) 
Praro = eee Patn = 


Observe tunnel temperature...T,,, = 


Ensure the following prior to tunnel runs: 
a. Top rack signal conditioners - ON (switch/light is on) 


Caution 
Do not change or attempt to adjust any settings at this time 


b. Traverse is in y = 0.00" position (visual inspection) 

c. Rotator is in 6 = 0° position, toggle switch = RUN, speed control = 75 
(visual inspection) 

d. 486DX/33MHz PC tumed ON and in Windows Program Manager 

e. HP Laser Jet II Plus turned ON with paper loaded 

f. Double doors and door to room H-033A CLOSED 


Double-click “Build-Time" under the "Wind Tunnel" menu (you are now in 
LabTech Notebook’s ICONview - you may need to click the uppermost right button 
to stretch the window so that you can see the entire set-up) 


Note: Because of the memory usage involved, waiting for WINDOWS 
to complete each input before assigning another 
will reduce the potential to overload the PC; if the PC 
is overloaded, chances are you will have to reset and start over. 
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5. Perform ZERO and SPAN 
Calibration 
The lab technician will ensure the data acquisition lab has the proper initial settings and 


that the input devices are properly calibrated. The configuration outputs should be 
"scaled" at this time using the set-up "ZERO" and "SPAN." 
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TABLE 3. AA2801 Cylinder Lab Initial" Settings 
a ee ee ee 


nea 
i 
5 


aa 

| | 
[iss>2| 19 | 
| Pac | m0 | 
a 

x100 

| eats" | Oem 
eee 


[Bir 


Bipolar 


‘Span () 


Sets 2 fe | 


Bipolar Bipolar 


"12 Bit Resolution" 
(2'? = 4096 divisions available) 
example: 


+1.25 - ( - 1.25) volts 
(4096 - 1) 


= + 0.305 mV resolution 
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“Table 3 show “initial settings" only and is intended for the technician to use prior to the 
beginning of the AA2801 Cylinder Lab. LabTech Notebook (LTNB) setups "ZERO," 
"SPAN" and accompanying directions should be used for final settings. 


Using "ZERO" and "SPAN" 


All input devices are referenced to some zero and to some secondary value, which yields 
a "span." That is, some mode of their output where the user wants 0.00 volts and some 
other voltage as a return. Everything is considered linear after the zero and span are set 
to a desired output. Normally a Digital Volt Meter (DVM) is used to verify such an 
output. However, for automatic data acquisition, it’s more important to verify the output 
at the PC since the date are being collected here. (It is possible for noise to affect what 
is read after the signal conditioner so that the signal has changed by the time it is 
displayed on the monitor.) The "ZERO" and "SPAN" set-ups were designed to give the 
user a visual display of his/her zero and span for each device used in the cylinder lab. 


Figure 22. ZERO Showing Signal Figure 23. ZERO Showing Signal After 
Somewhere Less Than Zero Volts Being “Flown” Up to Zero Volts 


Here are the steps: 


1. Using the 486 PC on the right, RECALL ZERO in LabTech Notebook Build-Time. 
Ensure the tunnel is not running by verifying the tunnel water micromanometer 
does in fact read 0.00 cm with the bottom of the water meniscus just resting on the 
cross-hair. Verify at the micromanometer, not the PC. 

3. Select RUN. 

4. The technician has already pre-configured the signal conditioners to yield a ballpark 
zero. Using the display presented, unlock and turn each "zero" knob on the upper 
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signal conditioners (#1, 2 and 10) until the colored "meatball" is centered. This 
is very easy for Naval aviators! A line display reflects the actual zero line for each 
device. Because of the complexity, the technician will ensure the rotator and 
traverse are properly scaled. 


5. To illustrate the effect of undesired noise on a data acquisition system, turn the 
DVM on and plug the red probe into the green outlet of the second signal 
conditioner; plug the black probe into the white outlet. Wait for 30 seconds. What 
happens to your display? (Don’t forget to unplug the DVM probes!) 

6. Exit ZERO and repeat step 1 for SPAN. 

7. After the operation brief and safety check of wind tunnel, run the tunnel up to 15 
cm H,O. 

8. Select RUN. This time, the signal will be referenced to 1.0 (unity). 

9. Again, the technician has already pre-configured the signal conditioners to give 
ballpark voltage output which, when scaled by the set-ups in LabTech Notebook, 
will yield 1.0. Using the display presented, unlock and tum each "span" knob on 
the upper signal conditioners (#1, 2 and 10) until the colored “meatball" is again 
centered. Another line display reflects the actual 1.0 line for each device. Again, 
the technician only will ensure the rotator and traverse are properly scaled. 

10. Press "REV" on the rotator control box to reset the cylinder to 6 = -20°. 

11. The setups are now properly scaled. Leave the tunnel running and proceed to step 
1 of the cylinder lab. 

Lab Procedures 

RUN #1 

1. RECALL PRE DIS1 

2. Single-click the "Run" menu and single-click "Run" 

(you will see a window which asks you if you want to "Overwrite previous 
log(s)?”) 

3.  Single-click "Yes" and expand the window (the display configured for this third of 
the lab is now present - data acquisition has been "triggered" to begin when the 
input from the azimuth block passes through -17° for the first time) 

4. | Momentarily push the black "FWD" button on the rotator control box, observe the 


following: 

a. _ the cylinder rotates from -20 to +180 degrees 

b. as the cylinder passes -17 degrees, data is displayed on the PC: 
Azimuth and current time displayed in meter boxes 


X-Y plot showin 
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Green line: theoretical C, vs. @ 
Cyan line: actual raw C, vs. @ (collected at Sampling Rate of 10 Hz) 
Black line: actual smoothed C, vs. 6 (moving average with parameter 
of 10 points) 
c. boundary layer separation occurs (the wake begins) 
d. data acquisition stops as cylinder passes through +180 degrees but the 
display remains on 


5. To make a print-out of the display hit ALT-Print Screen (this sends a BitMap to 
the ClipBoard in Windows’) 

6. Single-click “File” then single-click "Exit" to leave the display and return to 
ICONVview (you will be asked "Do you really want to exit?" - single-click "OK") 

7. Exit ICONview by single-clicking "File," "Exit," and "OK" in succession 

8. Double-click "paintbrush" in the Windows” Accessories menu 

9.  Single-click "View," then deselect (by single-clicking) "Tools and Linesize,” and 
"Palette" 

10. Click the uppermost right button to stretch the paintbrush window to its limits 

11. Select "Edit" and "Paste" (a copy of the previous display will appear) 

12. Select "File" and "Print" 

13. Reduce the scaling to 80% and hit “enter" (the printer will now print a copy of the 
display) 

14. Select the uppermost "File" and "Exit" (paintbrush will ask you if “you want to 
save current changes? - single-click "No”) 

RUN #2 

15. Repeat steps 1,2 substituting WAKE_ANI for the first set-up 

16. Single-click "Yes" and expand the window (the display configured for this third of 
the lab is now present - data acquisition has been “triggered" to begin when ALT-1 
is pressed on the keyboard...press ALT-1 when the next step is started) 

17. Have someone begin cranking the traverse handle at a quick, even, steady rate and 


observe the following: 

a. __ the traverse moves from 0 to 20 inches 

b. after ALT-1 is pressed, data begin to be displayed on the PC: 
Position and time are displayed in meter boxes 


X-Y plot showing 


Blue line: smoothed 


(moving average with parameter of 25 points and 
collected at Sampling Rate of 10 Hz) 


‘the position and relative strength of the wake 
d. data acquisition continues as traverse passes through 19 inches and display 


e 


remains on 

18. Press ALT-Print Screen to save display, then select File and Exit immediately to 
stop data acquisition. To make a print-out of the display repeat steps 5-14 

19. Stop the wind tunnel at this time using tunnel operating procedures, open test 
section 

RUN #3 

20. Repeat ZERO to recalibrate your signals, if necessary. 

21. Apply two 1/2" x 20" x 1 mil (approx.) pieces of cellophane tape to cylinder 1% 
inches on either side of the static port (this will trip the boundary layer, creating 
turbulence and thus delaying boundary layer separation) 

22. Crank the traverse back to y = 0 and close up the test section 

23. Run the tunnel up to 15 cm H,O again 

24. Repeat SPAN to recalibrate your signals, if necessary. 

25. Repeat steps 15-19 substituting WAKE_AN2 for the second setup (notice that the 
wake size is reduced) 

26. Note the tunnel temperature at shut down: T = 

27. You are complete - Exit LabTech Notebook and double-click "Exit Windows"; 

28. Place a formatted disk in either floppy drive. If you are using a 3'2" disk, at the 
DOS prompt type (otherwise substitute A for B): 

C:\> copy C:\windtunl\pre disl.m_B:\pre_dis1.m (ENTER) 

C:\> copy C:\windtunl\wake_ anil.m_ B:\wake_anl.m (ENTER) 

C:\> copy C:\windtunl\wake_an2.m_ B:\wake_an2.m (ENTER) 

This will copy the three data files onto your disk for later analysis. 
29. Secure the CPU, Monitor and Printer; ensure the wind tunnel is properly secured. 
Report Requirements 


1. 


Using a data analysis program such as MATLAB?’ and data from the lab, calculate 
the integrals of equations (1) and (2). Be sure to apply proper wake and solid 
blockage correction factors to find experimental drag coefficients for each of the 


three runs. See pages 35-38 of the NPS Laboratory Manual for Low Speed Wind 
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Tunnel Testing for more information. Explain in detail how you found the drag 
coefficients. List any equations used and how you used them. 


Using the c,, from Run #1, calculate the effective Reynolds number (Regrs) and 
determine the corresponding empirical drag coefficient using Figure 24 (from 


Anderson’s Fundamentals of Aerodynamics, p. 229) and 
Re rrr as Re ucorreceeds| +€) TF (8) 


Briefly discuss the differences between the pressure distribution results, the wake 
analysis results and the empirical drag coefficient. 


Determine the percentage reduction in drag afforded by an induced turbulent 


boundary layer. List at least two real-life examples of how turbulence is used to 
delay the onset of drag (boundary-layer separation). 
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FIGURE 3.39 
Variation of cylinder-drag coefficient with Reynolds number. (Experimental dasa as compiled in Panton, 
Ronald, Incompressibie Flow, Wiley-Interscience, New Yerk, 1984.) 


Figure 24. Cylinder Drag Coefficient 
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APPENDIX C: UPGRADED LAB DATA 


MMM MM MMMM ALAS AALS 


Re MME Re Re, 


Na tN MAE MMM MMU MUMMARBLARARAA A RAR ee RACARARAaeeaeaneeee.: 


OO TO ORI ON I RK 


WAKE ANI.m 


Filename 
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2 eA RRR eles. 
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APPENDIX D: MATLAB CODE FOR RETRIEVAL/DATA ANALYSIS 
Filename: DRAG1.m 


% DRAG1.M needs 3 files on a 3-1/2" disk labeled: 
B:\pre_dis!.m; B:\wake_anl.m; B:\wake_an2.m 
% If you haven’t placed a disk in the B-drive with these three files, as 
% named, do so at this time. \f you have our files ona 5-1/4" disk, 
u will need to get into this program * RAGI1.m" and edit the lines 
ith the load commands to read A vice B, for A-drive. 
emember, the student version of MATLAB will not execute a file with 
than 1024 elements. The Wake An files have at least that many. 
ext, it loads these three files into MATLAB and then begins to 
them down into column vectors for use in a modified Trapezoidal 
Instead of having equispaced x points, as required by the 
trapezoidal rule, the difference between each successive set of x-valucs 
are multiplied by the average between each associated set of y-valucs. 
These products are then summed. All cd’s are then corrected for blockage. 
The following abbreviations are used: 
cdpdu = = uncorrected cd ad we (for use > Re(cff) calculation) 
Reeff == effective Reynolds r for use w/ Anderson's Figure 3.39 
cdpdl == cd from pressure distribution method; eqn | basen serydoagent: 
cdpd2 == cd from pressure distribution method; eqn 2 - modi peat mthd. 
phil == 2 FOU f cree eR valucs was donc and a cqn | was 
cdwal == coefficient (RUN #2) via the wake analysis ans 2 
cdwa2 == ee coefficient (RUN #3) via the wake analysis method; eqn 2 
sepaz == the azimuth in degrees where separation occurs (estimate) 
hk re gmat alee dar eue an. 
== the wake (in y) 
redu == the % drag reduction afforded by ing the boundary layer 
sect some plots after the program finds <. ed. it any key to continuc 
after the computer Bene with each plot or cd. At this time type: 
drag! al oe 2 > oe Watch the B-drive light to 
clue you as to the proximity to your first value. It takes about 73 seconds. 
B:\pre_disi.m 
B:\wake_anl.m 
load B:\wake_an2.m 
pee = size(pre dist); 
{m2,n2}=size(wake_anl); 
{m3 ‘n3]= wee a2); 
x1 =pre_disl(1:m1,T); 
yl =pre_disi(1:m1 al 2); 
x2=wake_ani(1:m2,1); 
yeewes ani (1:m2 12); 
x3 =wake_an2(1:m3,1 a 
y3=wake an2(i:m3,2 
xx) =x1. *[(pi/180); % Neaiies to radian conversion 
y! ped ie -*cos(xxl); % cach elemental Cp is multiplied by its associated 
cos(thets 
% approx. equispace by letting delta x = total | # ray ie 
% Syprox. cigs by leg eqn. 5-16, p. 1 N; I call this 
edpdl = ((yy1 (m1)-yy1(1))/2 + (sum(yy1)- eb eet: yy: 
% this is oaly slightly diff fferent from eqn. I do not rely on 
% equispaced, x-values; instead | sum each cucesive product of dx and dy; 
% 1 call this “eqn. 2°: . 
for k=2:ml; 
80(K) = (yy LE *) + yy 1(k))/2*abs (xxl (k)-xx1(k-1)); 


en 
cdpd2 = sum(s0); 
p! ape ore) pi % this makes a polyfit of the values given to a Sth 


a 


ig et abet pg ea 


EEE 


‘ degree 

A emt % this evaluates that polynomial at each x-valuc 
p Al = (fib): Hives onthe (f(a) *(xxl(ml)/(ml-1)); % eqn 1 
or i=2 


2m2; 
sa 81 (i) = (y2(i-1) + y2(i))/2*abs(x2(i)-22(i-1)); % eqn 2 
cdwal = sum(s1)*2/2. 96; 
for j=2: 
324) = (36-1) + y3@)/2%abs(x3(j)-x3G-1)); % eqn 2 


ead 
cdwa2 =sum(s2)*2/2.96; 


cdpdu =cdpd2; 
cdpdu 
Lael =270477.9*1.049(1 + (.00234 + .0329*%cdpdu)); 
ce: 
% this next calculation adds the wake and solid blockage factors 
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nar =cdpdl *(1-3%((0.52°195.43)/(8.8 9144)" 1.5)-2°((2.96%cdpdl )/(2%45))); 
cdpd2 = cdpd2 *(1-39((0.52°195.43)/(8.89144)* 1. 5)-2° (2.96 %cdpd2)/(2%45))); 
c 

of = pftl *(1-39((0.52°195.43)/(8.8°144)* 1 .5)-2°((2.96 "phi )/(2°45))); 

P - 


plot(xxl yy] ,xx1,f) .grid,tithe(*Cp*coa(theta) vs theta’), xlabel(‘theta’),gtext(’Sth order polynomial’) 
cdwal =cdwal *(1-3%((0.52°195.43)/(8.89144)*1.5)-29((2.96*cdwal)/(2%45))); 
cdwal 
% what is the azimuth (location) where separation occurs? 
for t= 1:ml, 
if y1@+3) > y1@; 
sep = xx1(t);break ond 


end 
yp =2*2.96*sin(sep); 
sepaz =sep*180/pi; 


sepaz . 

plot(x1,y! ,sepaz,y1(t),"®’),grid,title(’Cp vs theta’), xlabel( theta’) ,gtext(’separation’) 
use 

wakebeg = 10-(yp/2); 

wak: 

wak =10+(yp/2); 

wakecad 


pause 
plot(x2,y2),grid,title("q ratios diff vs y*),xlabel(’y’) 

use 
cawad =cdwa2*(1-39((0.52%195.43)/(8.8%144)*1.5)-2*((2.96*cdwa2)/(2*45))); 
Cc 


percredu = 100-(cdwa2/cdwal)*100; 
percredu 


pause 
plot(x3,y3,x2,y2),grid,title("a comparison’), xlabel(’y’),gtext('Tripped Boundary Layer’) 
pause 
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APPENDIX E: MATLAB CODES FOR TUNNEL VELOCITY TABLES 


Filename: SPceEDS.m 


% This Matlab file SPEEDS.m creates conversion cohumas in different 
% wre bas on dat oui daring cabbratc the equations here 
% arc based on data obtained calibrations made by Professor 

- Hi ron 26 April 1993. 


This program is written oul to show conversions in ’ 
1 page formats and will pause and wait for any keystroke 
between displays to allow user to print data 


RRARRAA 
Py Fad 
PH a 
s 
at 
x a 
zs 
& 
t 
g 


RR 

oO 

m 
m 

7 

z 

8 

eG 


m/s kts) os mph = kav/hr 

is .05:1.0; % This is the ref.; cm of H20 on the micromanometer 

ne? ead: 9(62.35/30.5)-.24326))/(0. 89467 *.002378)); 

nor sqrt((2*(i. °(62.35/30.5)-.24326))/(0.89467 *.002378))*. 3048; 

mee sqrt((2°(i. 9(62.35/30.5)-.24326))/(0.89467*.002378))*.5921; 
rare Saqrt(24(i. %(62.35/30.5)-.24326))/(0.89467®.002378))*.6818, 

. un earn" %(62.35/30.5)-.24326))/(0.89467*.002378))°1.097; 

s= be f m’ &’ mi’ km’! 

Pi fi) grid xlabel("H20 cm’), ylabei('ft/sec") ,tuthe(" velocity’) 

Pet m) grid, xlabel(°H20 cm’), ylabel(‘m/s’),title( "velocity ’) 

Prat en. xlabel("H20 cm’), ylabel(’kts’),tithe('velocity") 

Plt. i xlabel(’H20 cm’), ylabel(‘mph’),title( velocity’) 


plot(i.km) grid, xlabel(” H20 cm’), ylabel(‘km/hr’),title( velocity") 
pause 


en 


(21.0.5: 10.0; % This is the ref.; cm of H2O on the micromanometer 
en f Faaqetaii *(62.35/30.5)-.24326))/(0.89467 *.002378)); 

“es ne 2 sqri(ant. (62.35/30.5)-.24326))/(0.89467 *.002378))*.3068; 
es k =sqrt((2°(i. %(62.35/30.5)-.24326))/(0.89467* 002378))*.5921, 


in 
m= sqrc(2(i. *(62.35/30.5)-.24326))/(0. 89467 *.002378))*.6818; 
% wn km/hr: 


ee ©(62.35/30.5)-.24326))/(0.89467 *.002378))*1 .097; 
oid f m’ k’ mi* km’} 


plot(.f), grid, alabel(*H20 cm’), ylabel('f/sec’),tithe(’velocity’) 
plot(ivm),grid,xlabel(H20 cm’), ylabel(’m/s’),titie(’velocity") 
plot(-k) grid, xlabel(*H20 om’), ylabel(’kts’),Wile(’velocity’) 

plot(ismi) grid, xlabel(H20 om’), ylabel(‘mph’),title(*velocity’) 


pause 

plot(i,km),grid,xlabel("H20 cm’), ylabel("ksxv/hr’),title(’velocity’) 
pause 

end 


I=100.: 3:20.0; % This is the ref.: cm of H2O on the micromanometer 
ee Feaqeuar6. . *(62.35/30.5)-.24326))/(0.89467*.002378)); 

_ m=aqrt((2%(i.*(62.35/30.5)-.24326))/(0.89467*.002378))*.3048; 
as East (62.35/30. 5)-.24326))/(0.89467*.002378))*.5921; 


94 


Baia (62.35/30. S)-.24326))/(0.89467°.002378))*. 6818; 


=agru2*4. v( *(62. 35/30.5)-.24326))/(O0. 89467 *.002378))* 1.097; 
peal te &’ mi’ km’) 


plott.f) grid, xlabel(H20 cm’), ylabel(’ft/sec’) ,wtke( velocity") 
plot(i a1) grid, xlabel("H2O cm’), ylabel(’m/s’), title(’ velocity’) 
plots £) grid, xlabel("H20 cm’), ylabel("kts’),tithe(‘velocity’) 
pause : . 
plot(i,mi),gnd,xlabel('H20 cm’), ylabel(’mph’),title(’ velocity ') 
pause 
plot(i,km), grid, xlabel(’H20 cm’),ylabel( km/hr’), tithe( ‘velocity ’) 
pause 
end 
® 
i= 20.0:.5:30.0; % This ts the ref.: cm of H20 on the micromanometer 
% in fUsec: 
f= sqrt((2°(i. *(62.35/30.5)-.24326))/(0.89467 *.002378)); 
% wm m/s: 
pole sqrt((2%(i. $(62.35/30. 5)-.24326))/(0.89467 *.002378))*. 3048; 
ein 
k=aqr((2%(i. %(62.35/30.5)-.24326))/(0.89467*.002378))*.5921; 
nin 
maar (62.35/30. 5)-.24326))/(0.89467*.002378))*.6818,; 
eum 
km = sqrt((2*(i. %(62.35/30.5)-.24326))/(0.89467*.002378))*1.097; 
Lo F m'k’ mi‘ km’) 
loti, f) grid, alabeK"H20 om'),ylabei(’ ft/sec’), tithe(’velocity*) 
peuse 
plot(i,m),gnid,xlabel("H2O cm’), ylabel(’m/s’), tithe(‘velocity’) 
pause 
plot(,k) gid, xlabel(’H20 cm’), ylabel(’kts'),tithe(’velocity’) 
pause 
plot(),mi),gnd,xlabel(°H20 cm’), ylabel(‘mph’),title("velocity’) 
pause 
plot(i,km),gnd,xlabel("H20 cm’), ylabel(’kin/br’),tidle(’velocity’) 


pause 
end 
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Filename: USER_SPD.m 


bs This Matlab file USER _SPD.m allows a user to provide MATLAB with 
% a desired velocity and fave MATLAB return a for the water 
% micromanometer (in cm H20) to be used with the kom Speed, as 


% calibrations made by Professor R. Howard and LCDR Clay: Maier oa 
% 26 April 1993. 


| AR 


roe MPH? y/n [o):’,’s’); 
if i= 
i = pati Eos desired with a decimal point: 
DEL_P=((k./.6818)*2%. eeer®: (002378/2)+ 94326). $162.35); 
("You need to set’), DEL 
FPo.s:.s 
meanad(d Ge. 35/30.S-0.24326))/(.895 *.002378))*.6818; 
plop. DEL P.k,’*’),grid,xlabel(’H20 cm’), ylabel(’ MPH’), title("velocity’) 


Pree’ Fusec? y/n (n]:’,’s”); 

if i= 
eae Enter speed desired with a decimal ) 
DEL_P=(k.*2*(.89467°*. eal ET "$162. 35); 
("You need to set’), DEL _P 


I=0: 5:.5:40.0; 
meqriteegee?. 35/30.5-0.24326))/(.895*.002378)); 
plod p.0 DEL _P,k,’*),grid,xlabel(H2O cm'),ylabel(’FUsec'),title(’velocity’) 


Break end 
elseif i= ='n’ 
ro m/s? y/n [n):’,"s’); 


if i rie = 
a desired with a 
DEL” veins " 5 8)*2*(.89467%. Oe 794326)°G0. 5162.35); 
(You need to set’), DEL_P 
J=0: 0.5:.$:40.0; 
aqrt((2*(j*62. 35/30.5-0.24326))/(.895 *.002378))*.3048; 
plo p.DEL_ P,k,’*’) grid, xlabel("H20 cm’), ylabel(’m/s’), title(’velocity’) 


ca end 
elseif i= ='n’ 
i= teat ak y/a {a]:’,’s"); 
if i ori =" 


"Enter speed des: 
DEL _ P=(k./.5921)*2%. Bac OOD 734326200. 5/62.35); 
(‘You need to set’), DEL_P 


j=0.5:.5 
p= grisea. 35/30.5-0.24326))/(.895 *.002378))*.5921; 
plotG,p,DEL_P,k,’*’), grid, xlabel(’H20 cm’), ylabel("kts’) stitle("velocity’) 


Preakad 


Lapel i=='n’ 

*kavVhr? y/n [n):’,’s") 

if i i=='n' Seats 
k=mput('Enter des wil 1 point: 
DEL_P=((k./1.097)*2%. oo 00237872) + "94326)°0. 5162.35); 
(‘You need to set’), DEL_P 


pause 

j=0.5:.5 

p= anita G62. 35/30.5-0.24326))/(.895*.002378))*1 .097; 
plotG,p,DEL_P,k,’*"),grid,xtabel(H20 cm’),ylabel(’kn/hr’), litle(’velocity’) 


ecorgey 
else 1= ='n' 
end 
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